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We analyze the semiclassical Hoi'ava-Lifshitz gravity for quantum scalar fields in 3 + f dimensions. 
The renormalizability of the theory requires that the action of the scalar field contains terms with 
six spatial derivatives of the field, i.e. in the UV, the classical action of the scalar field should 
preserve the anisotropic scaling symmetry {t — >■ L^^t, x — > L^x, with 2; = 3) of the gravitational 
action. We discuss the renormalization procedure based on adiabatic subtraction and dimensional 
regularization in the weak field approximation. We verify that the divergent terms in the adiabatic 
expansion of the expectation value of the energy-momentum tensor of the scalar field contain up to 
six spatial derivatives, but do not contain more than two time derivatives. We compute explicitly 
the counterterms needed for the renormalization of the theory up to second adiabatic order and 
evaluate the associated /3 functions in the minimal subtraction scheme. 
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One year and a half ago, Hofava proposed a new approach to formulate a quantum theory of gravity pj]. Hofava's 
theory, which has attracted enormous attention, consists of a non-diffeomorphism-invariant ultraviolet modification 
" of Einstein's general relativity. The main idea in [l[ is to extend Einstein- Hilbert action with higher spatial derivative 
^ ' terms, whose introduction, while manifestly breaking local Lorentz invariance, leads to heal the short distance diver- 
^ ^ I gences and ultimately yields a power counting renormalizable theory. The way this is achieved without introducing 
ghost instabilities is keeping the requirement of the theory to be of second-order in time derivatives. This introduces 
' an asymmetry between the time coordinate t and the coordinates associated to a preferable foliation that de- 
' , fines a three-dimensional space-like hypersurface of induced metric '■'^'sij. In turn, four-dimensional diffeomorphism 
invariance results manifestly broken at short distances, and consequently the theory only exhibits diffeomorphism 
invariance in three-dimensions, in addition to the reparameterization invariance in time. According to this picture, 
"*7^ ■ the four-dimcnsional general covariance of gravity would emerge merely as an approximate symmetry at low energy. 

Fragmentation of space-time diffeomorphism invariance in the form of a preferable three-dimensional space-like 
^ • hypersurfaces defined at constant time, immediately suggests to consider the ADM decomposition for the metric as 
the convenient picture. Namely, consider 



^ , ds^ ^ -N'^df + g,j{dx' + N'dt){dx^ + N^dt), (1) 

>: 

' where, as usual, Latin indices refer to the spatial coordinates, i,j — 1,2,3, and gij = ^^^gij- In the non-projectable 
C theory, the lapse function N depends both on time and the spatial coordinates, in such a way general relativiy is 

5^ captured within this formulation. 

The action of Hoi^ava's theory is given by 



O 



S ■ 



[ dt dx^ Ny^ {K^jK'^ - XK^ -2A + ^R~V) (2) 



O . 167rG 

where A is the bare cosmological constant and A and ^ are arbitrary bare coupling constants; Einstein the- 
.'^ ory corresponds to the special choice A = 1. The extrinsic curvature Kij in the ADM variables takes the form 
Kij — + \/iNj + \/ jNi)/ {2N), whose trace is given hy K = Kijg^^ . Here, Vi denotes the spatial derivative, while 
^ ■ the dot denotes the derivative with respect to time. The function in © plays the role of a potential, as it only 
depends on spatial derivatives of the metric, which would include higher derivative contributions (see below). In ([2]), 
R represents the Ricci scalar curvature of the three-dimensional space-like hypersurface of induced metric gij . 

The presence of terms in the action that involve higher spatial derivatives leads to different scaling dimensions for 
the time and the spatial coordinates. This is represented by the scaling symmetry 



t L''- t. 



N,. 



(3) 
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which is characterized by the dynamical critical exponent z. A consistent choice is z = 3, which is the one we will 
consider throughout this paper. With this choice, we can consider the potential 

V = i(aiAi? + a2R^JR'^ + ...)+ 47rG(6i A^i? + b2R,jR^''Rl + ... ), (4) 

where we are using the notation A — ViV\ The ellipses in ^ stand for other terms of the same dimension. 

At low energy, the action turns out to be dominated by the term that involves the Ricci scalar R, with coefficient 
^. In turn, the theory would reproduce Einstein's general relativity in the infrared, provided A flows to the value 
Air = 1. The consistency of the theory and the validity of this hypothesis were extensively discussed in the literature; 
see 0-13 and references therein. Of special interest is the discussion in [1, 0] , where an improved version of Hof ava 
gravity, which seems to be free of pathologies, was presented. 

About renormalizability, of particular importance is the question on how the coupling of Hofava gravity to matter 
affects the properties of the theory in the UV. With the purpose of addressing this problem, we study the coupling of the 
theory to a quantum scalar field, representing the matter content. The gravitational field will be treated at a classical 
level, so we are considering a semiclassical Hofava-Lifshitz gravity. It is interesting to remark that, if the matter fields 
satisfy the usual dispersion relations (i.e if the classical action has four-dimensional general covariance), the theory 
is non rcnormalizable. Indeed, it is well known in the context of quantum field theory in curved spacetimes that in 
order to absorb the divergences associated to the matter fields it is necessary to include in the gravitational action 
terms proportional to TZ^ ,TZ^^TZ'^'^ and TZ^^palZ'^"'"^ , where TZ^^pa denotes the components of the four-dimensional 
Riemann tensor. These terms contain four time derivatives of the metric, and therefore are not included in Hofava 
gravity. As we will see, renormalizability of the field theory demands that the action for the matter sector contains 
terms with six spatial derivatives, implying that in the UV the coupling to the scalar field preserves the Lifshitz-type 
anisotropic scaling with critical exponent z — 3. We will verify that the divergent terms in the adiabatic expansion of 
the expectation value of the stress-tensor associated to the scalar field actually contains up to six spatial derivatives 
but it remains of second order in time derivatives. We will explicitly compute the counterterms needed for the 
renormalization of the theory up to second adiabatic order, and we will write down the corresponding /3-functions in 
the minimal substraction scheme. 

The computation techniques we will employ here have recently been employed with success to study renormalization 
in the so-called Einstein-aether theory and in other field theories with modified dispersion relations Q. The idea for 
using the same techniques in Hofava-Lifshitz gravity comes from the observation that this theory is closely related to 
such Lorentz violating scenarios; see for instance Q and [loj . 

Let us begin by considering the coupling of Hofava-Lifshitz gravity to a Lifshitz-type scalar field. In the ADM 
form, the components g^^i, of the four-dimensional metric ([1]) are given by 

where i,j = 1,2,3, and gij refers to the metric on the three-dimensional foliation of constant t. We consider small 
perturbations of the metric about fiat space; namely, we write 

N = l + dn, N'^6N\ g,j^d^j+hij. (5) 

We consider a matter Lagrangian giving by a scalar field ip that also exhibits anisotropic critical scaling; namelv(TTj| 

where the potential F{ip,dip) is given by 

F{ip, dip) = -gid''p}d^Lp - 52 {^^f + 93^'^'P^f, 

where we have to be reminded of the definition Aip — -^di{y/gg^^ djip), with g^^ = — {'^)g^i — ((^)(;*" ^^-'g-'")/^''''!?'^'^, 

and (^)(7°'' = —N^^. The equation for the Green function reads 

-dt [^{dt - d,)G{x,x')^ + dj {^^^{dt - N'^ d,)G{x,x')^ - m^^N G{x,x') + 2g^d, {N^ d'G{x,x')) - 

252V5A(iV /\G{x, x')) + gs^A {N A^G{x, x')) + g^^A^ {N AG{x, x')) = ~5{x - x'). (6) 

At the linearized level we have ([5]), which yields ^ =1-1- ft./2, g^^ — 6ij — hij. This can be used to write the 
equation for the Green function ([B]) in the weak field approximation. The Feynman propagator of zero order in the 
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metric perturbations reads 



G^p^ {x, x') — 
where /c = |/c| and 



(2^)4 (^-kl+ul^zs)' 



ul^m^ + 2gik^ + 2g2k'^ + 2g3k^ , 
while the first order contribution can be written as follows 



d'^k f d^p e^P-^e''^-^^-^') /fc(p) (1 + ep)"^ 



G^p\x, x') 



{2n)^J (27r)4 (^k^ + - ^e) {-k^, + u^l) ' 
Here, fkip) is a function of kQ,ki,pQ, and that is linear in the metric perturbations 

fk{p) = (Sn -'^)kl + 2ikihSN' - (Sii + ^)ujI + hjk'k^'^ + h,jSrsk'k' k'-p' 



(7) 



+ik4Pa \^6n ^ J ^ SN^poki + ik^SN^pi — ( ( + 77 ) % p'fc^ + h^^Pikj 



(5 71 



d(fc2)3 \^ 4 
d3.;2 



p'k' + lh,,k^k^{Srsk^pn' - ^-^{k.k^fk^ 1^ (^/k,kY - \h^,ky 
^^p'^k^S.jkY + {S.jkYf - \hij¥yp^5rsk''p' ~ \hijp'k^{5rsk''p'')' 



1^ + ^P' (%fcVr - ^/^.-fc^fc^/ - lh.AsPYk^k^P' 

|^(^P%fcV-^/..fcV/). (8) 



with p2 = |pj2^ g^j^j^ (;jefine(;j a.s 



= .2^ '2 ^ ■ 9 



In what follows, for the sake of convenience, we perform a Wick rotation in Eq.© and we call k^ = ikg. 

To obtain the adiabatic expansion of the Feynman propagator we start by expanding the integrand of Eq.Q in 
powers of po and pi. The different adiabatic orders of fk{p) sue given by 



ff^"^ = {Sn-'-^)kl + 2ikikiSN' - (5n+^)ujl + hijk'y^, (10a 



duji 



dk^' 



+h,,5rsk'yk-f-^^, (10b) 
(t^'^' + ^/».,FF(<5.,fcV)2 - ^(5,^.A:V)2A:2) . (10c) 
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From Eq.® one can see that no powers of po appear in the adiabatic orders /^'^'''"^ with to > 3. Besides this property, 

these m > 3 adiabatic orders are not relevant for our discussion. Note also that the adiabatic orders f^'^'""^^ with 
TO > 6 vanish. 

We have also to expand ep in its adiabatic orders. That is, 

,ad(i) _ 2^fc4Po + 2^ %fcy 

P ^ ujl + kl ' ^ ' 

" ujl + kl ' ^ ' 

and = 0. It is easy to see that ep'^'-™-' with to > 3 do not involve powers of po! explicit expressions of these 

adiabatic orders are not necessary for our present purposes. 

The expression from which one obtains the components (T^i/(a;)) of the expectation value of stress tensor after 
taking the coincidence limit x a;' is obtained by evaluating the derivatives of the propagator Gf{x,x'), as usual. 
Appropriate regularization turns out to be necessary. For the case of the component Tqq (x) , we find 

(Too(x)) = \mildtdt. + -{SN'dA'+SN'd,,dt) + —{l + 26n) 

- (1 + 25n) (^-g,6''' d,dl - g^d^d'^ + \9z{d^d'^ + d^d'^) 

+ ^53 {hijdidjd'^ + hi,j>dedj>d^ + diKjdjd"^ + dih^j'dj'd^) 

- 93 (h''^' d^d^,dj,d''^ + h'^d'^d.djd'^ ~ dji^jd^d'^d''^ - d^h^j,d^dj,d''^ 

- ^93 (d^h'^d.djd'^ + d^h''^' d^d,>df - d.d^h.^djd'^ - d,d%j>d^dj. 

- 93 (dkh'J d^d.dkd'^ + dk'h''^'d^d,,dj,dk' - dAhjdkdjd'^ - d,dk-hrd''dk'd,,) } Im Gf{x, x') (12) 

where hij = hij — ^Sij, = didi^ and a primed index on a derivative indicates that the derivative is taken with 
respect to a primed coordinate. 

For the sake of simplicity, and because it is enough for our present purposes, we partially fix the gauge by setting 
(5iV* = 0. This greatly simplifies the expression of the (Toi(a;)) component, which reads 

{TM) = \ lim {dtd^, + dt'd^) Im Gf(x,x'). (13) 

Z x—^x' 

To obtain the regularized expectation values of the stress tensor we use dimensional regularization. Therefore, 
after computing the derivatives of hnGp^x^x') that appear in Eans. (|12p and (jl3p . we can set x = x' . Then, 
it is straightforward to separate the different adiabatic orders of {To^{x)) (/i — 0,1,2,3), before performing the 
integrations. In this way we obtain an integral expression for each adiabatic order. The next step is to use dimensional 
regularization to perform the integrals. We apply dimensional regularization to both temporal and spatial directions, 
but in a separated way. That is, we split the d— dimensional integrals into integrals in di and d2 dimensions, with 
d = di + c?2, where di — >■ 1 and d2 3. All integrals in are of the form 

where j is an even number {Id^{j,l) = if j is odd) and / is an integer; = /T{di/2), with T{z) being the 

Gamma function. Note that the right hand side is finite in the limit c?i — > 1 (as is usual in dimensional regularization 
for odd dimensions), then, we set di = \. 

Using the adiabatic expansion of the expectation value of the stress-tensor, by simple power counting one can study 
up to which adiabatic order it contains divergences and how many temporal derivatives do appear in the divergent 
terms. To illustrate this, let us consider as an example the following contribution to {Tqq): 

twix) = lim 9t9t'Im Gp{x,x'). 
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It can be shown that this term, along with others which are similar, are the most divergent ones. The contribution 
that is linear in the metric perturbations can be written as 

(27r)4" J i2nr {kl+ulf 



Let us first analyze the terms that do not contain pq. In a schematic way, it is simple to show that the ultraviolet 
behavior of a term in fkip) and in Cp, respectively, is given by Cp ^ k^~"p" /{k^ + w^), with 1 < n < 6, and 
fk{p) ~ 5gk^~'^p'^ ^ where < s < 5 and 5g represents a component of S^'^'^g^^. Then, for a term characterized by n, s 
and r, the integration in k^ yields 



{kl+Lolf^ ^' +W^)2+'- V/c/ ^ \k) 

Therefore, by power counting it can be shown that the integral in k is convergent only if s + rn > d + 2. That is, for 
d = 4, we have that — is the maximum power of p that appears in a divergent contribution; i. e., terms of 
adiabatic order greater than six are finite. 

Let us now analyze the terms that contain powers of po- In Ea. (jl4p po appears explicitly and also implicitly through 

fkip) and ep. Notice that of aU the adiabatic orders only /^'^'"^^(p), ep^^^^ and ep''^^'' depend on po (see Eans. lfTOl) and 
dTTI) and the paragraphs that follows each of them) . Then, one can write all the terms of second adiabatic order that 
involve one or two powers of po (i-e., terms with poPi or with p§ and no additional power of and show that all of 
them are logarithmically divergent. Hence, as the convergence improves with the adiabatic order, we can conclude 
that the contribution of higher adiabatic orders with at least one power of po will be finite. Below, we compute 
explicitly the terms of second adiabatic order of the 00 and Oi components of {T^y) and we show that while the terms 
with one and two powers of po are both logarithmically divergent, the ones with Pq do not appear in the final result. 
We expect that terms with p^ do appear in the ij components, but here we will not compute these explicitly. 

In summary, in {T^^) there appear divergences up to in the sixth weighted adiabatic order, where the weighted 
adiabatic order of a term is given by (z = 3) 

W ~ zriQ + Tii 

where Ui and uq are, respectively, the number of spatial derivatives and time derivatives appearing in the term. This 
is analogous to the weighted power counting criterion introduced in [l2l| for field theories in Minkowski spacetime (see 
also [13]). 

In order to illustrate the procedure by which we obtain the regularized adiabatic orders in terms of 3— tensorial 
quantities, let us consider as an example the zeroth adiabatic order of {Tqq(x)). After performing the integrals in pp 
(which are straightforward) and the integral in k^ (as described above), we obtain 

{TooixT"'^"^ = 4~ / (0^ {} + ^'^'^ " ^) " ^ + 252fc.fc,fc2 + 3g,hk,k') | , (15) 

where fj, is an arbitrary parameter of dimensions of mass, introduced to ensure that ip has the correct dimensionality. 
To carry out the angular integrations we use the following property [l^ : 



d'^^kk'K..k''-g{k^) = 



if r is odd, 

T'^-'-Arig] if r is even. 



where 



T'l-^'- = hs'^'^5">'^...S''-^''- +a\\ permutations of the i's]. 



r 



oo 



r[l/2]r[(d2 + r)/2] J„ 



The remaining integrals can be related by performing an integration by parts, 

,4— do /■+00 



,,4— dr) r+oo //4— dr) 

w" 
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where we have defined Iq = dk'^k'^'^ ^w^. Moreover, one can easily show that, due to the gauge condition 

5N'^ = 0, one has {Toi{x))'^'^^^^ = 0. Therefore, as expected, the lowest adiabatic order of the energy momentum 
tensor is proportional to the metric and can be absorbed into a redefinition of the cosmological constant (see below) . 

We follow the same procedure for the second adiabatic order of (To/i {x))- After a long but straightforward calculation 
the results are 



^2Go. + , ^, {2Go^ + Zm]) \ , (18) 



where Goo and Goi are components of the linearized Einstein tensor G^^ and i^j is the linearized extrinsic curvature, 
and we have defined the following integrals: 







t^fe dk? ' 7o ' Jo ^1 d{k'^Y 



Notice that while Ii is quartically divergent, I2 and are logarithmically divergent. The terms in Eq. p^ have first 
time derivatives of the metric and, as we have anticipated, result to be logarithmically divergent. We have repeated 
all the calculations without partially fixing the gauge SN"^ — 0, and reobtained Eq. ([T8| as a cross-check. 

With these results we can now analyze the renormalization of the bare constants associated to the terms of second 
adiabatic order that appear in the gravitational action To do so, we start by writing the 00 and Oi parts of the 
semiclassical equations for the metric (in the weak field approximation), keeping only terms up to second adiabatic 
order; namely 



StJg {^^'^30^ + ^^™} = = (^oo(x)).e„ + (Too(x))^'iW + (roo(x))-'i(2)^ (^9) 

^ {Go. - (A - l)d,K'^ = (To,(x)) = (To,(x)),e„ + (To,(x))'^<i(2)^ (20) 

where we have added and subtracted the adiabatic expansion of {To^{x)) in order to separate the renormalized part 
{Toi_i{x))ren and the divergent contributions. The latter are to be absorbed into a redefinition of A, G, A and ^. Then, 
we introduce Eqns. (HZ]) and ^ into (HH) and and we find that (Too(x))'''^(°) and (To^(x))'''*(2) can be cancelled 
with the following choice of the bare constants: 



AG- = (AG-),-^i^-P^/o, 



G-i 



6(27r)'*-2 



6(27r) 



I2- 



2h 



(d~l)(d+l) 



G-i(A-l) = (G-i(A-l))fl,+ 



2(27r)'i-2 (d- l)(d + l)' 



(21a) 
(21b) 
(21c) 
(21d) 



where we denote the renormalized constants by a subscript R. 

It is worth noting that from these equations we can recover the well-known results corresponding to the usual 
{z — \) scalar field by setting gi — 1/2 and 92 = 53 = before taking the limit d — )■ 4. In such a case, vanishes 
and in the limit d — > 4 we have Iq ^ m^{d — 4)^^/4 and /i = /2 ~ m? /{d — 4). 

All the integrals on the right hand side in Eqs. ((2T|) are divergent in the limit d ^ A. In the particular case of a 
massless field (m = 0) these integrals can be computed explicitly. We assume that 33 > and — 43153 > in order 
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to avoid zeros of ujk. Thus, in the limit d — >■ 4 we have: 



l2 



^3 



h = - 



{2n) 

(27r)<^-i 
^i^-^^d-i 



(2tt 



■ -u, — ± J 

h 
h 



92(9^ - 4.91.93) 
{9I - 4gig3) 



1 . , 1/4 



d-A 
1 



1 



1 



2 

572^ 



rf-4 
1 



d-4 

111(^53^") 



ln(A'53^'') 



FP, 



d-A 



where FP denotes the /x-independent finite part. 

The renormahzation group equations are obtained simply recalling that the bare constants are independent of /i, 
and are given by 



92 {9I - 4ffig3) 

1 (g|-4.9ig3) 

\/2 127rg^/' 
3 



5v2g^7r 



(23a) 

(23b) 

(23c) 
(23d) 



Note that, in contrast to what happens in the case of a usual (z = 1) scalar field, here we have obtained that the 
renormalized constants depend on /i for a massless field. Note also that in these equations, as f/3 > and g2—'^9i9d, > 0, 
the right hand side have a determined sign, except in the first equation which depends on the sign of 52 • 

Assuming that the order of magnitude of the constants gt are determined by a single mass scale Mq, i.e. gi ~ 
q^Mn ~ gsM^ ~ C(l), these equations give us information on the running of the renormalized constants in the UV 
[16| . at scales larger than Mc- With the renormahzation scheme we are considerin g (m inimal subtraction), it is not 
possible to analyze the IR behaviour of the coupling constants. As pointed out in |14|. [TtI this would require a 
mass dependent renormahzation scheme. We hope to address this problem in a forthcoming publication. 
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